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Abstract We discuss a model for the autocatalytic reaction A + B ^ 2 A 
on substrates where the reactants perform a compact exploration of the 
space, i.e., on lattices whose spectral dimension d is < 2. For finite systems, 
the total time r for the reaction to end scales according to two dilTerent 
regimes, for high and low concentrations of reactants. The functional de- 
pendence of T on the volume of the substrate and the concentration of reac- 
tants is discussed within a mean-field approximation. Possible applications 
are discussed. 



1 Introduction 

Diffusion-reaction processes is a long standing problem which finds a number 
of applications, especially in physics [IJ, chemistry and biology 2 . 

Most of the earlier studies focus the attention on a single particle diffus- 
ing in the presence of immobile reactants, while much less is known about 
the statistical properties associated with the diffusion of a set of particles, 
notwithstanding its interest. Indeed, multiparticle diffusion problems are 
difficult to manage due to the fact that the effects of each single particle do 
not combine linearly, even in the noninteracting case [5]. In the last years 
much effort has been devoted to the formulation of rigorous many-body 
treatments of diffusion controlled reactions especially in low dimension. In 
fact, while in high dimensions a mean-field approach provides a good de- 
scription, in low dimension local fluctuations are responsible for significant 
deviation from mean-field predictions [1] . 

In general, a great deal of recent experimental as well as theoretical 
work has been devoted to the study of such diffusion-reaction processes in 
restricted geometries. The latter expression refers to two different (possibly 
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concurrent) situations: that of a low dimensionality, and that of a small 
spatial scale. In the first case, the spectral dimension d characterizing the 
diffusive behavior of the reactants on the substrate is low (1 < c? < 2), and 
the substrate underlying the diffusion-reaction lacks spatial homogeneity. 
Hence, there are considered media whose properties are not translationally 
invariant and where the reactants perform a "compact exploration" of the 
space [5] . These kinds of structures can lead to a chemical behavior signif- 
icantly different from those occurring on substrates displaying an homoge- 
neous spatial arrangement. This is the case, for example, of fractal lattices: 
in the last 20 years, an extensive literature has been investigating the conse- 
quences of a fractal geometry on the laws of reaction kinetics |6] , for example 
for the one-species A ^ 0)[7] and two-species {A + B ^ 0) [7(l8l [9] anni- 
hilation reactions. In all these systems the role of the generally noninteger 
spectral dimension, whose definition will be discussed below, is stressed, 
as opposed to the integer euclidean dimension characterizing homogeneous 
structures. 

But restricted geometry also refers to a variety of experimental situa- 
tions in which these processes occur on spatial scales too small to allow an 
infinite volume treatment. The so-called finite-size corrections to the asymp- 
totic (infinite-volume) behavior in this case become predominant. Indeed, 
previous works considered infinite systems (both euclidean and fractal), 
and studied their properties in some kind of thermodynamic limit; typi- 
cally, sending the volume to infinity while keeping the density of reactants 
fixed. Therefore, they considered the critical properties of the systems (for 
example, the scaling of the density of reactants) for t oo, hence, for an 
infinite time lapse of the reaction. One of the most important issues of this 
paper regards the finite size of the systems under study. In this work we 
examine explicitly finite systems where no thermodynamic limit has to be 
taken. All the quantities we calculate, in particular the total reaction time 
r, are hence finite, and we seek their dependence on the finite parameters 
of the system (volume of the reaction and concentration of the reactants) . 

In particular, in this paper we study the dynamics of a system made 
up of two species particles undergoing irreversible quadratic autocatalytic 
reactions according to the following scheme: A + B ^ 2A, with reaction 
probability set equal to one. All particles move randomly and particles of 
different kinds react on encounter, i.e. the reaction is strictly local and de- 
terministic. Autocatalytic reactions have been extensively analyzed on Eu- 
clidean structures, both analytically and numerically PII1[TT1[T^131I141I15] . 
A continuous picture of this system can be attained by the Fisher equation 
[16pi7j which describes the system in terms of front propagation; however, 
this picture will not intervene in our calculations, that will mostly concern 
the low-density regime, where a front propagation cannot be defined. 

While previous works on autocatalytic reactions considered reactions on 
Euclidean lattices, here we focus, as mentioned above, on low dimensional 
structures (1 < d < 2), hence considering media whose properties are not 
translationally invariant. Our investigations are especially concerned with 



Autocatalytic reaction on low-dimensional substrates 



3 



the role of topology in the temporal evolution of the system. In particular, 
we will consider the concentration pA{t) of A particles present in the sys- 
tem at time t and its fluctuations. From pA{t) it is also possible to derive 
an estimate for the reaction velocity. Furthermore, we consider the average 
time T at which the system achieves its inert state, i.e. Na = N. We call 
this time "Final Time" . As we will show, r depends on the number of par- 
ticles N and on the volume V of the underlying structure, meant as the 
total number of sites. More precisely, it will be shown, both numerically 
and analytically, that for small concentrations of the reactants the "Final 
Time" factorizes into two terms depending on N and V respectively. This 
results agrees with previous works where the model under study was ana- 
lyzed for Euclidean lattices [T8lll9j. We will also show how this dependence 
could provide a practical tool for the determination of the concentration 
of reactants, especially when very small concentrations have to be detected 

m- 

The plan of the paper is the following. In Sec. [2] we introduce the model, 
we recall previous results on Euclidean substrates (Sec. l2.l( ) and discuss the 
main features of inhomogeneous lattices (Sec. 12. 2|) . In Sec. Owe show our 
analytical results concerning those lattices; Sec. [4] discusses the results of 
numerical simulations. Sec. [H] contains our conclusions. 

2 The model 

We consider a system made up of N particles of two different chemical 
species A and B, diffusing and reacting on a discrete substrate with no ex- 
cluded volume effects. The volume of the substrate is V; at time t, NA{t) and 
Nsit) are the number of A and B particles, respectively, with N = Na+Nb- 
We define pa(0 = NA{t)/V and psit) = NB{t)/V as the concentrations of 
the two species at time t. 

Different species particles residing at time step t, on the same node or 
on nearest-neighbor nodes react according to the following mechanism: 

A + B ^2A 

with reaction probability set equal to one, so that the process is strictly 
diffusion-controlled. Notice that the previous scheme is quite general as it 
also includes possible additional products (other than 2A) made up of some 
inert species of no consequences to the overall kinetics. 

The initial condition at time t = is NaIO) = 1, Nb{0) = N~1, with aU 
particles distributed randomly throughout the substrate. As a consequence 
of the chemical reaction defined above, NA{t) is a monotonic function of t 
and, due to the finiteness of the system, it finally reaches value N; at that 
stage the system is chemically inert. The average time at which NA{t) = N 
is called "Final Time" and denoted by r. 

The Final Time r is of great experimental importance since it provides 
an estimate of the time when reaction-induced effects (such as side-reactions 
or photoemission) vanish [21] . 
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In this perspective, deviations from the theoretical prediction of r are, 
as well, noteworthy: they could reveal the existence of competitive reactions 
[22] or explain how the process is affected by external radiation [23] , 

However, one of the most interesting applications of the Final Time is 
analytic [20ll24j : as we will show, r sensitively depends on N, that is on the 
initial amount of reactant. Hence, given a trace reactant, its determination 
can be achieved by means of spectrophotometric measures of t. 

Indeed, our results confirm that this technique can be extremely sensitive 
pnil^ and the determination of ultratrace amounts is therefore possible. 

We stress again that he following analysis mainly focuses on high-diluted 
finite-size systems. The experimental applications previously described con- 
cerns finite systems; furthermore, we aim to evidence the role of the sub- 
strate topology which just emerges in the diffusion-limited regime. 



2.1 Euclidean Substrate 



The quadratic autocatalytic system for diffusing reactants on Euclidean 
lattices has been analyzed in detail in earlier works [TSlfTO] in the context of 
information spreading. It also provides a simple model for epidemic systems: 
A (B) particles stand for (irreversibly) sick (healthy) or informed (unaware) 
individuals respectively. For these systems as well a knowledge of the rate of 
infection or information diffusion is of great importance. We briefly review 
the results obtained in p^lT9] for r on Euclidean lattices. 

In general, r depends on system parameters N and L and, in the low- 
concentration regime, and for the Euclidean lattices this dependence can be 
factorized into two contributions depending on N and L respectively, and 
whose functional form depends on the dimension of the lattice. A mean-field 
calculation for r provides exact results for d > 2 and c? = 1 , while dimension 
d = 2 is marginal. Our previous results can be summarized as follows: 



t{N,L) 



r2 



d = 1 



f2{N)LHii{L) d 



(1) 



2 + HN)^d 

N 



L'^ d>3, 



where Cd are dimension-depending constants and fitting procedures suggest 
/2(iV) = ii+Mi^ (for further details see ll9j). 



2.2 Inhomogeneous Structures 

A number of problems in physics and chemistry are related to random walks 
on homogeneous structures. In the last few years there is a growing interest 
in the theory of diffusion in low and fractal dimensions [26]. In fact, most of 
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the materials existing in nature show a disordered, non-crystalline geometri- 
cal structure. Indeed, fractals are good models for disordered systems since 
they usually display a dilation symmetry which is a fractal characteristic. 
Apart from the applications point of view, low dimensional systems are of 
great theoretical importance. As previously underlined, when diffusion is 
considered, low dimensional structures dramatically affect the kinetic laws. 

Fractals are self-similar structures exhibiting dilation symmetry. Dif- 
ferently from Euclidean structures, their description requires at least two 
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Fig. 3 Comb Lattice 

(typically noninteger) different dimensions. The first is the fractal dimen- 
sion df, that gives the dependence of the volume of the system (i.e., the 
number of sites it contains) on its linear size 

V{L) - L'^f. 

The second is the spectral dimension d, which governs (among other phe- 
nomena) the long-time properties of diffusion on the lattice. Indeed, if we 
consider a random walkers starting from a given site i of the lattice, the 
probability Pii{t) of returning back to the starting point at time t, at large 
times, follows the law 

Also, the number S{t) of different sites visited by the random walker at time 
t is for large times: 

S{t) - t^/"^ for d < 2 
S{t) - t for d > 2 (2) 

For d < 2 the random walker returns to its starting site with probability 
1 and the lattice is called recurrent] for d > 2 the probability of return is 
less than 1 and the lattice is called transient (lattices with d = 2 have to 
be discussed case by case). For d < 2 the walker is also said to perform 
a compact exploration of the space [5] , since the (fractal) dimension of the 
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random walk trajectory is greater than the dimension df of the underlying 
lattice. 

The spectral dimension is of a more general interest than the fractal 
dimension, since it can also be defined for inhomogeneous structures that 
lack a dilation invariance, hence for which a fractal dimension df cannot be 
defined. 

The two fractals we consider in this paper, the Sierpinski gasket and the 
T-fractal, have both d < 2. 

The Sierpinski gasket is generated by iterating in a recurrent fashion a 
generating cell consisting of a triangle (Fig[T]). The number of iterations is 
called the generation g of the fractal. The total number of triangles after g 
iterations is 3^, while the total number of sites (hereafter called volume V) 
is V = 3^-^. The linear size of the gasket is given by 2^ . The Sierpinski 
gasket has fractal dimension df = « 1.584 and spectral dimension 

J= i2g9 « 1.365. 

logs 

The T-fractal is constructed from a 4-sites T-shaped generating cell 
(FiglH). It has fractal dimension df = « 1.584 and spectral dimen- 

sionJ=iH||«i.226. 

The third structure we will consider, the comb lattice (Figin]), does not 
present dilation invariance; hence, the df cannot be defined. As we said 
above, it is still possible to define a spectral dimension that turns out to be 
J= 3/2. 

Some linear problems have already been solved exactly on these struc- 
tures using renormalization groups technique (see for example 27J). How- 
ever, for many interacting diffusing particles an exact solution is not feasible 
and we rely mainly on numerical simulations. 



3 Analytical Results 

In this section we study the irreversible autocatalytic reaction occurring in a 
close system by means of a mean field approximation. This kind of approach 
is very different from those previously adopted for this kind of system. 

An analytical result for the dependence of the Final Time r on V and 
N is difficult to find: approximate calculations can be carried out in the two 
limit regimes of high and low concentration. 

For high concentration (p ^ 1) the results found in Euclidean structures 
[18] continue to hold for inhomogeneous geometries: the A particles occupy 
a connected region of the space for all t. The border of this region expands 
at a fixed velocity, such that at time t the region covers all the sites whose 
chemical distance from the starting point of the seed particle is <2t. Hence, 
for a finite system like those we are considering here, the Final Time is 
T = ^max/2, where Zmax is the chemical distance of the most distant point 
on the lattice, starting from the seed particle. On Euclidean geometries 
this yields r = L/2ford=l and r = L for d > 2. On the other hand, 
on inhomogeneous structures the dependence on L is not so simple, since 
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it involves taking the average with respect to aU possible starting points 
for the seed particle; anyway, the crossover between this high-concentration 
regime and the low-concentration one remains apparent. 

Our mean-field approach is based on a different point of view and the 
assumptions introduced make it valid just in the low-concentration regime. 
In this approach we focus on collective quantities lacking the spatial depen- 
dence. 

In particular our hypothesis is that the time elapsing between a reac- 
tion and the successive one is long enough that the spatial distribution of 
reactants can be considered as uniform. This assumption corresponds to a 
mean field approach since we neglect correlations between spatial position 
of reactants; in other words we neglect the fact that for small times after 
a reaction the two A particles are likely to be find nearby. This kind of 
hypothesis is therefore valid for small concentration of reactants, i.e. p <C 1. 
As a consequence of this hypothesis we can just focus on two-body interac- 
tions among particles since the event of three or more reactants interacting 
together is unlikely at small densities. 

First of all we consider the final time t. Let us define (i„) the time 
elapsing between the n — 1-th first encounter among different particles and 
the n-th one. This time corresponds to the average time during which there 
are just NA{t) = n particles in the systems. In the mean field approximation 
this is proportional to the trapping time in the presence of n traps randomly 
distributed through a volume V. For compact exploration of the space {d < 
2) , the average trapping time ttrap for a random walker in a distribution 
of iV — Na randomly distributed moving traps is given by [28j 




since the density of traps is ptraps = NA{t)/V. Here, the symbol ~ denotes 
proportionality. 

This is the trapping time for one particle in a background of moving 
traps; we are interested in the average trapping time of the first out of 
N — Na particles, that for rare events is just the same time rescaled by a 
factor N — Na (the number of B particles) : 





{t„) = V^^'^^ 



2/d 



N - Na 



(3) 



The time t can therefore be written as 



N-l 




(4) 



Na = 1 



For this sum there exists no closed form; however, in the limit N ^ oo we 
can adopt a continuous approximation and replace the sum with an integral. 
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Fig. 4 (Color online) Concentration of A particles pA{t) vs time t for a system 
made up of = 128 particles embedded on different structures, as shown by the 
legend. The best fit for the cubic lattice is also shown and it is given by a pure 
sigmoidal. 



We then find 



V 



2/d 



(2 - d)N 
where is the harmonic number 

Hm — 



(5) 



fc=l 



In particular, the leading-order contribution for a one-dimensional system 
(J = d = 1) is 

(6) 



N 

and for a two-dimensional lattice (c? = d = 2) 

JogA^ + 7 



V- 



N 



(7) 



where 7 is the Euler-Mascheroni constant. For 1 < d < 2, the expression in 
^ interpolates between © and ([7]). 

The mean-field extension just performed also allows to derive some in- 
sights into the temporal behavior displayed by NA{t). In fact, being t„ the 
average time at which Na = n, from Eq. [4] we can write 



Tn= ^^^^^ = 

Na = 1 



(8) 
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Now, we estimate NA{t) as 

whose numerical solution provides an S-shaped curve to be compared (see 
Fig. IH) with the sigmoidal curve obtained from a standard mean-field ap- 
proximation [18j 

A'' 

^A(t) = (^_,)^-.,.^, . (9) 

where p is a quantity proportional to the concentration p that in practice 
must be adjusted within the fitting procedure. 



4 Simulations 

In this section we show results obtained with numerical simulations per- 
formed on the Sierpinski gasket, on the T graph and on comb structures. 

First of all we consider the dependence on t displayed by NA{t). In Fig. [5] 
and[5]we show data obtained for the Sierpinski gasket. In particular, in the 
latter figure we also provide a comparison with results obtained for the T- 
graph, the comb lattice, the square lattice and the cubic lattice. Consistently 
with results found in [TS], for transient lattices N^it) is well fitted by the 
sigmoidal function of Eq. [HI while for low-dimensional structures deviations 
are expected. 

In fig|4]we show the time evolution NA{t) for substrates with the same 
total number of particles N and with (approximately) the same volume V, 
but different spectral dimension d. In order of decreasing d, they are: the 
cubic lattice {d — d — 3), the square lattice {d = d = 2), the comb graph 
(J = 3/2), the Sierpinski gasket fractal (^ = i§|f « 1-365) and the T- 
fractal {d — |^|-| « 1.226). We remind that the spectral dimension describes 
the long-range connectivity structure of the substrate and the long-time 
diffusive behavior of a random walker on the substrate. In particular, from 
Eq. [21 we expect that for substrates d < 2 the number of different sites 
visited by each walker will grow faster as d increases, and so will the number 
of meetings between walkers. Hence, we expect the curves NA(t) to grow 
faster, and saturate earlier, with increasing d {N and V being fixed). This is 
precisely what happens, as shown in Fig. [4| (except for the saturation time 
T of the comb lattice, which will be discussed below). For d > 2 (e.g., d — 3 
in the figure), NA{t) is independent of d and is fitted by a pure sigmoidal 
function. 

From NA{t) one can derive the rate of reaction 

v{t)^dtNA{t). (10) 

As you can see from the numerical results in Fig. [6l in agreement with the 
theoretical one, v(t) is an asymmetrical curve exhibiting a maximum at 
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Fig. 5 (Color online) Data depicted in the figures above refers to a system of 
A'^ = 128 particles diffusing on a Sierpinski gasket of five different generations. 
Each generations is depicted in different colors, as shown by the legend. Left: 
Number of A particles NA{t) present in the system vs time t. Right: Fluctuations 
a{t) — {N\) — (Na)"^ as a function of time 
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Fig. 6 (Color online) Reaction velocity v{t) for a system of N = 128 particles 
diffusing on a Sierpinski gasket; five different generations (depicted in different 
colors) are shown. The reaction velocity is defined in Eq. [TD] 



a time denoted by t^. This time obviously corresponds to a flex in N^lt) 
which scales with the volume of the structure according to the following: 

This is the same dependence shown by r (see below), and corresponds to 
a situation in which the population of the two species are about the same 
{Na = Nb = N/2). Analogous results can be obtained for the T-fractal. 

Furthermore, the profile shown in Fig. [6] also suggests that the efficiency 
of the autocatalytic reaction is not constant in time but, provided the num- 
ber N of particles is conserved, it exhibits a maximum when the number of 
B particles is about N/2. 
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Fig. 7 Scaling of r with the linear size of the system for a 1-dimensionaI chain 
(blue circles), a Sierpinski gasket (black triangles), a T- fractal (red squares), and 
a comb lattice (green plus signs) on a double-logarithmic scale. The number of 
reactants is fixed at A'' = 1024 for all systems. Dotted lines highlight the low- 
concentration regime (L 3> 1), corresponding to a power law for all systems. 
For the 1-dimensional chain, the linear high-concentration regime (low L) is also 
highlighted. 



A similar result may be derived for the variance a a (t) of the number of 
A particles present on the substrate. 

Interestingly, fluctuations display a maximum at a time which, again, 
depends on the system size with the same law as r. Notice that t^r > ty and 
NAiU) = f . 

Finally, we consider the time r representing the average time at which 
the autocatalytic reaction stops since all B particles have been transformed 
into A particles. In general this quantity depends on system parameters V 
and N = Na + Nb and, as we will show, its functional form is significantly 
affected by the topology of the substrate. 

The dependence of r on the system size L (or the volume V — L'^f ) 
clearly displays two different regimes, as shown in Fig. [71 In both cases, r 
increases with L, but in the high-concentration regime the growth is less 
rapid; in particular (as shown for d = 1 in the figure), it is proportional to 
L for Euclidean lattices. 

In the low-concentration regime, where we can assume that reactions 
only occur among two particles, the mean-field- like calculation explained in 
the previous section holds and we expect r to vary with N and V according 
to Eq. [5] which can be rewritten as: 



(11) 
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Hence in Figure [8] we plotted tV '^^'^ vs N and we fitted data according 
with the r.h.s . of the previous equation. 

It can be seen that for small densities all the data collapse. Moreover, 
in that region, the fit coefficients introduced are in good agreement with 
theoretical predictions. 

It should be underlined that the average final time depends non-trivially, 
on N and V, viz. t does not depend directly on the total concentration p, 
though the dependence on N and V can be factorized. 

The agreement of formula [H] for the comb lattice is less good. In partic- 
ular, it seems that the dependence of r on and V can still be factorized, 
but that the exponent for V is rather 8/3 than 2/c? = 4/3. This may mean 
that the particular mean field approximation we have made does not hold 
anymore for strongly inhomogeneous structures such as combs; this point is 
still under investigation. 

As explained in Section [21 experimental measures of r are useful in 
monitoring trace reactants. Indeed, our results show that r — fj^{N)g^{V) 
and therefore, once the substrate size is fixed, the initial amount of rcactant 
can be expressed as iV = /j"^( g-("v) )■ 

A proper estimate of the sensitivity of this method is provided by the 
derivative the smaller the derivative and the larger the sensitivity. In 
Fig. [H] we depicted numerical results for both N and its derivative ^ as 
a function of r; topologically different substrates are also compared, all 
sharing, approximately, the same volume. The numerical plots provided 
allows a qualitative analysis and comparison among the different structures 
considered; a more quantitative inspection can be outlined after a proper 
calibration procedure. 

First of all, notice that the characteristic curves are well defined and 
they allows a univocal determination of N from t. Moreover, the sensitivity 
of this analytic technique is better for small values of N (with N > 2) and, 
interestingly, for low-dimcnsional substrates. Indeed, once V is fixed, when 
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Fig. 8 (Color online) Rescaled Final time tV~'^^'' vs number of particles A'^ 
for Sierpinski gasket (top) and T-fractal (bottom). Different symbols and colors 
distinguish different generation as explained by the legend. The line provides the 
best fit in agreement with equation 1111 
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d < 2, the technique sensitivity is improved by lowering the substrate dimen- 
sion. On the other hand, when d > 2, the resuhing curves are overlapped, 
hence no improvement is achieved. It should be underlined that the high 
sensitivity attained just for small concentrations of reactants makes this an- 
alytic technique very suitable for the determination of ultratrace amounts 
of reactants, which is of great experimental importance [29t.30u31j 

5 Conclusions 

We have presented a model that considers the autocatalytic reaction A + 
B 2 A on non-Euclidean, low-dimensional {d < 2), finite-size substrates, 
characterized by a volume V and a total number of reacting particles N . 

We showed by analytical calculations that the Final Time r (the total 
time span of the reaction) displays two different regimes, for high and low 
concentrations, with a different dependence on V and N. In particular, the 
functional law for low concentrations can be recovered by means of a mean- 
field approximation. In fact, with respect to the standard one, our mean- field 
approach is able to take into account the topological effect arising from a 
low-dimensional substrate. 

Numerical simulations corroborated these results for fractals, while sim- 
ulations on strongly inhomogcneous lattices (combs) hint at a quantitatively 
different behavior. 

Theoretical results concerning the average Final Time find important ap- 
plications in analytical fields, where measures of r are exploited for detecting 
trace reactants. Our results suggest that the sensitivity of such technique 
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Fig. 9 Loglog scale plot of the reactant amount TV vs Final time r; as shown in 
the legend, different substrate structures (with approximately the same volume) 
have been compared. In the inset the derivative ^ is depicted again as a function 
of r. Lines are guides for the eyes. 
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is affected not only by the reactant concentration, but also by the topology 
of the structure underlying diffusion. More precisely, a small concentration 
of reactants implies a better sensitivity, hence allowing the determination 
of ultra-trace reactants. Moreover, at the same concentration, and for low- 
dimensional {d < 2) substrates, by reducing the dimension d the sensitivity 
can be further improved. 
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